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We explore the electronic properties of finite-length graphene nanoribbons as well as graphene 
nanodisks with various sizes and shapes in quest of metallic ones. For this purpose it is sufficient 
to search zero-energy states. We find that there exist no zero-energy states in finite-length zigzag 
nanoribbons though all infinite-length zigzag nanoribbons have zero-energy states. The occurrence 
of zero-energy states is surprisingly rare. Among typical nanodisks, only trigonal zigzag nanodisks 
have degenerate zero-energy states and show metallic ferromagnetism, where the degeneracy can be 
controlled arbitrarily by designing the size. A remarkable property is that the relaxation time is 
quite large in spite of its small size in trigonal zigzag nanodisks. 



I. INTRODUCTION 

Grapheneiii^ii^, a single atomic layer of graphite, has 
invoked a great interest in the electronic properties 
of graphene-related materials. In particular, graphene 
nanoribbona^i^i^i^'^'iSiiii'i^ have attracted much atten- 
tion due to a rich variety of band gaps, from met- 
als to wide-gap semiconductors. They are manufac- 
tured by patterning based on nanoelectronic lithogra- 
phy methodsi^iiiii^. It is interesting that graphene with 
a zigzag edge has the half-filled fiat band at the zero- 
energy level and show edge ferromagnetism^. The half- 
filled zero-energy states emerge also in all zigzag nanorib- 
bons and hence they are metallio^i^. However, realistic 
nanoribbons have finite length. It is important to investi- 
gate the finite-length effects on the electronic properties 
of nanoribbons. 

Another basic element of graphene derivatives is a 
graphene nanodiskl-. It is a nanometer-scale disk-like 
material which has a closed edge. A graphene nan- 
odisk can be constructed by connecting several benzenes. 
There are many type of nanodisks, where typical exam- 
ples are displayed in Fig[TJ Finite-length nanoribbons 
may be regarded as nanodisks provided that their length 
is short [Figl2]. Some of nanodisks have already been 
manufactured by soft-landing mass spectrometrj^i^. 

In this paper we analyze the electric properties of nan- 
odisks as well as finite-length nanoribbons. Since all 
zigzag nanoribbons are metallic, as we have mentioned, 
we expect all zigzag graphene derivatives are also metal- 
lic. On the contrary, the emergence of zero-energy states 
is quite rare. We show that there are no zero-energy 
states in finite-length zigzag nanoribbons. We also in- 
vestigate a class of trigonal and hexagonal nanodisks 
possessing zigzag or armchair edges, among which we 
have found zero-energy states only in trigonal zigzag nan- 
odisks. 

Trigonal zigzag nanodisks are prominent in their elec- 
tronic properties because there exist half-filled zero- 
energy states and they are metallic. Indeed, we can en- 
gineer nanodisks equipped with an arbitrary number of 
degenerate zero-energy states. Furthermore, spins are 
argued to make a ferromagnetic coupling. A remarkable 
property is that the relaxation time is quite large in spite 




FIG. 1: Basic configurations of typical graphene nanodisks. 
(a) Benzene, (b) Trigonal zigzag nanodisk (phenalene). (c) 
Trigonal armchair nanodisk (triphenylene). (d) Hexagonal 
zigzag nanodisk (coronene). (e) Hexagonal armchair nanodisk 
(hexa benzocoronene)^^. 

of its small size. 

This paper is organized as follows. In Section [ll] and 
\ni\ based on the tight-binding Hamiltonian, we calcu- 
late the energy spectra of finite-length zigzag nanorib- 
bons and of a wide class of graphene nanodisks, respec- 
tively. In Section we also carry out a systematic 
investigation of the zero energy states in trigonal zigzag 
nanodisks. In Section ITVl we analyze the wave functions 
of these zero-energy states to examine how they are lo- 
calized at the edges. In Section [Vl we study the spin- 
relaxation time of nanodisks. Section [Vl] is devoted to 
discussions. 



II. ENERGY SPECTRUM OF FINITE-LENGTH 
NANORIBBONS 

We calculate the energy spectra of graphene deriva- 
tives based on the nearest-neighbor tight-binding model, 
which has been successfully applied to the studies of car- 
bon nanotubesi^ and nanoribbons^. The Hamiltonian is 
defined by 

H = Y.^^4c^ + J2*^l''^^''l^ (2.1) 

where is the site energy, tij is the transfer energy, and 
cj is the creation operator of the tt electron at the site i. 
The summation is taken over all nearest neighboring sites 
Owing to their homogeneous geometrical configu- 
ration, we may take constant values for these energies, 
Ei = £f and tij = t. Then, the diagonal term in (|2.ip 
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FIG. 2: Geometric configuration of zigzag nanoribbons with 
width W and length L. Here we show the example of the 
{W, L) = (3, 8) nanoribbon. The basic chain is W connected 
benzene depicted in gray. Short nanoribbons may be regarded 
as parallelogrammic nanodisks. 
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yields just a constant, £f^c, where A^c is the number 
of carbon atoms in the system. The Hamiltonian (|2.ip 
yields the Dirac electrons for graphenei^^ii^. There exists 
one electron per one carbon and the band-filling factor 
is 1/2. It is customary to choose the zero-energy level 
of the tight-binding Hamiltonian (|2.ip at this point so 
that the energy spectrum is symmetric between the pos- 
itive and negative energy states. Therefore, the system 
is metaUic provided that there exists zero-energy states 
in the spectrum. 

In this section we investigate finite-length nanoribbons 
to know if there are zero-energy states. A classification of 
infinite-length nanoribbons is given in a previous work^. 
Here we concentrate on finite-length zigzag nanoribbons. 
We classify them as follows (Figl2]). First we take a ba- 
sic chain of W connected carbon hexagons, as depicted 
in dark gray. Second we translate this chain. Repeating 
this translation L times we construct a nanoribbon in- 
dexed by a set of two integers {W, L). In what follows we 
analyze a class of finite-length nanoribbons generated in 
this way. Parameters W and L specify the width and 
the length of nanoribbons, respectively. The infinite- 
length nanoribbons are obtained by letting L — > cx). The 
finite-length nanoribbons are regarded as parallelogram- 
mic nanodisks when L ^W. 

In analyzing a nanoribbon containing Nq carbon 
atoms, the Hamiltonian l|2.ip is reduced to an iVc x Nc 
matrix. It is possible to diagonalize the Hamiltonian ex- 
actly to determine the energy spectrum Ei together with 
its degeneracy gi for each finite-length nanoribbon. The 
density of state is given by 



(2.2) 



i=l 



It can be shown that the determinant associated with the 
Hamiltonian 112.11) has a factor such that 



det [si - H (Nc)] oc (e - t)<^'^)(£ + (2.3) 



implying the a (W, L)-fold degeneracy of the states with 



FIG. 3: (Color online) Density of states of finite-length 
nanoribbons. The vertical axes is the energy e in units of 
t — 3eV, and the horizontal axes is the degeneracy. There 
exist no zero-energy states. 



the energy e = ±t, where 



a{l,L) = 2,1,2,1,2,1,2,1,2,1,... 

a(2,L) = 1,1,0,2,0,1,1,0,2,0,... 

a(3,L) = 2,0,2,0,2,0,2,0,2,0,... 

a(4,L) = 1,2,0,3,0,2,1,1,2,0,... 



(2.4a) 
(2.4b) 
(2.4c) 
(2.4d) 



We have displayed the full spectra for some examples of 
finite-length nanoribbons by taking t = 3eV in FiglSl 

One of our main results is that there are no zero-energy 
states in finite-length nanoribbons. However, the band 
gap decreases inversely to the length, and zero-energy 
states emerge as L oo, as shown in FigHl This 
is consistent with the fact that infinite-length nanorib- 
bons have the fiat band made of degenerated zero-energy 
states^i^. Hence, a sufficiently long nanoribbon can be 
regarded practically as a metal. 



III. ENERGY SPECTRUM OF NANODISKS 

We next derive the energy spectrum of each nanodisk 
[Figll] by diagonaHzing the Hamiltonian (|2.ip . As an 
example we display the density of state l|2.2p of trigonal 
zigzag nanodisks in FigUl We have classified them by 
the size parameter N as defined in FigEJa) . The number 
of carbons are given by Nc — N'^ + QN + 6. 

In order to reveal a global structure, it is convenient to 
introduce the doped electron number at a given energy 
E, which we normalize as 



\E)^^f D{e)de, (3.1) 
^^c Jo 
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FIG. 4: Band gap of zigzag nanoribbons as a function of 
length L. The horizontal axes is the length L and the vertical 
axes is the energy gap AE in units of t = 3eV. Each curve is 
for width VK = 1 to 7 from top to bottom. 
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FIG. 5: (Color online) (a) Geometric configuration of trigo- 
nal zigzag nanodisks. It is convenient to introduce the size 
parameter A'' in this way. The 0-trigonal nanodisk consists of 
a single Benzene, and so on. The number of carbon atoms are 
related as A^c = N'^ + QN + 6. See explicit examples given 
in Figlll (b) Density of states of the A^'-trigonal nanodisk for 
= 0, 1, 2, • • ■ ,7. The horizontal axis is the size A'^ and the 
vertical axis is the energy e(A'') in units of t = 3eV. Dots on 
colored bar indicate the degeneracy of energy levels. 



with < 1. We then make the inversion of this 

formula to derive E as a function of n. See Figl6] for 
their correspondence in graphene, where E (n) is found 
to be a prolonged S-shaped curve. 

Diagonalizing the Hamiltonian (|2.ip explicitly, we have 
constructed and displayed E{n) for several nanodisks 
with trigonal zigzag shape in Figlllja), trigonal armchair 
shape in Figl7l[b) and hexagonal zigzag shape in Figl7l[c). 
Each diagram consists of step-like segments reflecting the 
i5-function type density of states l|2.2p . The length of each 
step represents the degeneracy of the energy level in units 
of Nq ■ It is remarkable that there exist zero-energy states 
only in trigonal zigzag nanodisks. We have also checked 
explicitly the absence of the zero-energy state in a series 
of nanodisks with hexagonal armchair type [Fig[lje)]. 

In each figure we have also displayed the prolonged S- 
shaped curve of graphene, which the E (n) of nanodisk 
approaches in the large size limit (A^c ^ oo). The pro- 
longed S-shaped curve is universal regardless of the nan- 
odisk's shape. 

We investigate trigonal zigzag nanodisks more in de- 
tails since they have zero-energy states. It can be shown 
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FIG. 6: (a) Density of state D(e) as a function of the energy 
e in units of t = 3eV for graphene. (b) The energy E{n) as a 
function of the doped electron number n. These two functions 
share the same information of the energy spectrum. 

that the determinant associated with the Hamiltonian 
(|2.ip has a factor such that 

det[e/-i7(iVc)] cxe^(e-i)°^^^(e + i)"^^\ (3.2) 

implying the iV-fold degeneracy of the zero-energy states 
and the a (iV)-fold degeneracy of the states with the en- 
ergy e = ±i, where 

a {N) = 3, 3, 3, 3, 3, 5, 3, 5, 3, 7, 3, 7, 3, 9, 3, 9, 3, (3.3) 

for iV = 1,2,3,- ■•. 

Since there exist half-filled zero-energy states for A^ > 
1, these nanodisks are metaUic. The gap energy between 
the first excitation state and the ground state decreases 
as the size becomes larger. However, it is remarkable 
that the gap energy is quite large and is of the order 
of 3eV even in the nanodisk with size N = 7 [FigE^b)]. 
This is much higher than room temperature. Hence, the 
low-energy physics near the Fermi energy £ = can well 
be described only by taking the zero-energy states into 
account. 



IV. WAVE FUNCTIONS OF TRIGONAL 
NANODISKS 

A wave function is represented as 

i 

where ipi{x) is the Wannier function localized at the lat- 
tice point i. The operator Ci in the Hamiltonian (|2.ip an- 
nihilates an electron in the state described by the Wan- 
nier function ipi{x). We are able to calculate the am- 
plitude uji for zero-energy states in the trigonal zigzag 
nanodisk. All of them are found to be real. As an ex- 
ample we show them with size A^ = 5 in FiglHl where 
the soHd (open) circles denote the amplitude uJi are pos- 
itive (negative). The amplitude is proportional to the 
radius of circle. It is intriguing that one of the wave 
functions is entirely locaHzed on edge sites for nanodisks 
with A^ =odd, as in FiglHJa). There are no such wave 
functions for nanodisks with TV =even. 
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FIG. 7: (Color online) Energy spectrum E{n) as a function 
of the doped electron number n for (a) trigonal zigzag nan- 
odisks, (b) trigonal armchair nanodisks and (c) hexagonal 
zigzag nanodisks, with various sizes. The horizontal axis is 
the number of doped electrons n, and the vertical axis is the 
energy E in units oi t = 3eV. The S-shaped solid curve is 
that of graphene. The energy density diagrams of nanodisks 
approach that of graphene for large size. There are degener- 
ate zero-energy states in all trigonal nanodisks, and they are 
metallic. There are no zero-energy states in all zigzag arm- 
chair nanodisks and all hexagonal nanodisks, and they are 
semiconducting. 



The most important property is that all wave functions 
are nonvanishing on edge sites. In order to demonstrate 
this, we have investigated how the zero-energy states are 
modified by changing the site-energy Si in the Hamilto- 
nian (|2.ip only for edge carbons. Edge carbons are those 
surrounded by two carbon atoms and one hydrogen atom, 
while bulk carbons are those surrounded by three car- 
bon atomsS.. If a wave function vanishes on edges, the 
zero-energy state must remain as it is. First, we take 
Ei = £ — Ae for all edge carbons. We show how the 




FIG. 8: The zero-energy states of the trigonal nanodisk with 
size A'' = 5. There are 5 degenerate states, (a) ~ (e). Elec- 
trons are localized on edges in the state (a). When the site 
energy Si is decreased at edges equally, the degeneracy is par- 
tially resolved, as illustrated in FiglOja). The state (a) has 
the lowest energy and nondegenerate; the states (b) and (c) 
are degenerate; the states (d) and (e) are degenerate and have 
the highest energy. When the site energy Si is decreased fur- 
ther on the bottom edge, all the degeneracy among these 5 
states is resolved, as illustrated in Figl9][b). 
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FIG. 9: (Color online) Zero-energy states split into several 
nonzero-energy states by edge modifications. The horizontal 
axis is the size A'^ and the vertical axis is the energy in units 
oit — 3eV. (a) We take Ei = e — Ae for all edge carbons with 
Ae = O.OSeV. Nonzero-energy states are nondegenerate (blue) 
or 2-fold degenerate (red), (b) We take Zi = e — Ae — Ae', 
Ae' — O.OOSeV for edge carbons on only one of the three edges 
and ei = e ~ Ae for those on the other two edges. All states 
become nondegenerate. 



zero-energy states split by taking Ae = O.Olt = O.OSeV 
in FigEJa) for TV = 1, 2, • • • , 20. It is seen that all 
zero-energy states acquire negative energy and that they 
become one-fold or two-fold degenerate. The energy 
decrease is larger when the wave function is locaHzed 
more on edges. When we decrease the site-energy more, 
ei = e - Ae - Ae', Ae' = O.OOlt = O.OOSeV for edge car- 
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bons on only one of the three edges, all the degeneracy 
is found to be resolved as in Figl9l[b). 



V. MAGNETIC PROPERTIES OF TRIGONAL 
NANODISKS 

We proceed to show that the degenerate ground states 
lead to a ferromagnetism as in graphene nanoribbon. 
This is because the Coulomb exchange interaction or the 
Hund's rule coupling drives all spins to polarize into a 
single direction. We are most interested how large is the 
relaxation time for small nanodisks. (Strictly speaking, 
ferromagnetism can occur only in an infinitely large sys- 
tem, but we may use the terminology for a finite system 
if the relaxation time is large enough.) 

The effective Hamiltonian for the Coulomb exchange 
interaction is given by the Heisenberg model. Here, for 
the sake of simpHcity, we use the Ising model. 



H 



N 



(5.1) 



since the Heisenberg model presents essentially the same 
result on the relaxation time as we shall argue later. In 
the effective Hamiltonian, ai is the spin operators of elec- 
trons in the ith zero-energy state, ai = ±1, and the 
summation is taken over all electron pairs. An impor- 
tant point is that the exchange interaction strength Jy 
must be nonzero for all electron pairs because their wave 
functions are nonvanishing on edges. This is in a sharp 
contrast to the Hamiltonian for a nanomagnet, where Jij 
can be regarded nonvanishing only for neighboring elec- 
tron pairs since the index i represents the site in the real 
space. This makes a clear difference in the relaxation 
time as we shall soon see. 

The spin dynamics is well described by the master 
equation, dP {t) /dt = £P{t), where the symbol P de- 
notes the probability distribution function specifying the 
spin configuration, and C is the Liouville operator as- 
sociated with the Hamiltonian. When the expansion 
P (0) — J2x ca-Pa holds at the initial state at time t — 0, 
the state of the system at any later time is given by 



P(t) -e^*P(0)-^CAe-^*PA 



(5.2) 



where A is obtained by solving the eigenvalue equation, 
£P\ = — APa- The relaxation rate of the system is equal 
to the minimum eigenvalue Ai„in of the Liouville operator. 
Thus the relaxation time is given by r = 1/Amin- 

To get a concrete idea, are nonvanishing, 

we first make an estimation by making an approximation 
Jij — J. We are able to diagonalize the Ising Hamiltonian 
(|5.ip explicitly. The eigenvalues are given by 
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FIG. 10: Relaxation time of various graphene nanodisks. 
The relaxation time of nanodisks with size A'' = 1, 2, • ■ ■ , 8 
from bottom to top. The horizontal axis is the coupling con- 
stant J/kT and the vertical axis is the relaxation time in the 
form of logio[Tjv/Ti]. (b) The relaxation time of nanodisks 
with interaction strength K = J/kT = 1/16, 1/8, 1/4, 1/2, 1, 2 
from bottom to top. The horizontal axis is the size of nan- 
odisks A''. 



where n = 0, 1, 2, • • • ,N is the energy level index. We 
have then carried out an exact diagonalization of the 
eigenvalue problem of the Liouville operator C and deter- 
mined the eigenvalue Amin- We show the relaxation time 
as a function of the coupling strength J/kT for various 
size N in FigfTOTa) and as a function of size for various 
coupling strength in FigfTOlfb). 

For the noninteracting case, J/kT — 0, it is easy to 
see that the relaxation time r^r of the iV-trigonal zigzag 
nanodisk is given by 

Tat — Nti. 

On the other hand, the relaxation rate Amin is given by 
the Arrhenius-type formula for strong coupling limit or in 
low temperature limit, J/kT 3> 1. Hence the relaxation 
time is given by 



TN — exp 



AE' 



Tl, 



where 



AE = JN^/2 for N = even, 

AE = J{N^ - 1) /2 for N = odd 



(5.4) 



(5.5) 



is the energy difference between the highest energy state 
and the ground state. It is observed in FigdOlthat the 
relaxation time is given by the asymptotic formula (|5.4p 
already for J > kT. 

The relaxation rate is given by the Arrhenius-type for- 



mula also in the generic model (|5.ip with Jij ^ J, where 
AE is the energy difference between the highest energy 
state and the ground state. Then, the relaxation time is 
given by (|5.4p by replacing J with 



^efF ■ 



1 



N 



N{N - 1) 



(5.6) 



which is of the order of a typical Jij . Furthermore, we 
would obtain the same result for the Heisenberg model 



(5.3) 

even with Jij ^ J. 
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Note that in an ordinary nanomagnet composed of N 
spins, the relaxation time is given by l|5.4p with 

AE = JNz, (5.7) 

where z is the number of the nearest neighboring spins. It 
is remarkable that the size dependence of the relaxation 
time is oc iV^ for trigonal zigzag nanodisks, though it is 
oc N for normal nanomagnets. This is because any one 
spin couples with all other spins in the zero-energy state. 

VI. DISCUSSIONS 

Graphene derivatives become metallic when they have 
half-filled zero-energy states. We have explored the en- 
ergy spectra in a wide class of nanodisks as well as finite- 
length nanoribbons. It is surprising that the emergence 
of zero-energy states is quite rare. There exist no zero- 
energy states in finite- length zigzag nanoribbons. How- 
ever, the band gap decreases inversely to the length, and 
zero-energy states emerge as L ^ oo. Hence, a suffi- 
ciently long nanoribbon can be regarded practically as a 
metal. 

Among a wide class of nanodisks we have studied, trig- 
onal zigzag nanodisks are distinguished for their elec- 



tronic properties since they exhibit metallic ferromag- 
netism due to their half-filled degenerate zero-energy 
states. The degeneracy is controllable arbitrarily by 
changing the size of nanodisks. We have estimated the re- 
laxation time, which has been argued to be proportional 
to exp[JeffiV^/2fcT] when the size is N. Though the nu- 
merical estimation of the effective spin stiffness Jeff is yet 
to be done, it is of the order of the Coulomb energy since 
its origin is the exchange interaction or the Hund's cou- 
pling rule. We conclude that the relaxation time is quite 
large for its size at low temperature T < Jeff/2k. Hence, 
for instance, it could be used as a memory device. By 
connecting nanodisks with nanoribbons, we can design 
electronic circuits. These devices would be obtained by 
etching a single graphene. Alternatively, nanodisks may 
be connected with leads by making tunneling junctions. 
We would like to make an analysis of the nanodisk and 
lead system together with related phenomena in future 
works. 
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